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Abstract The ACES (Atomic Clock Ensemble in Space) mission is an ESA - CNES project 
with the aim of setting up onboard the International Space Station (ISS) several highly sta- 
ble atomic clocks with a microwave communication link (MWL). The specifications of the 
MWL are to perform ground to space time and frequency comparisons with a stability of 
0.3 ps at one ISS pass and 7 ps at one day. The ACES mission has applications in several 
domains such as fundamental physics, metrology or geodesy. 

The raw measurements of the ACES MWL need to be related to the scientific prod- 
ucts including all corrections (relativity, atmosphere, internal delays or phase ambiguities) 
and considering all terms greater than 0. 1 ps when maximized. In fact, the mission aims at 
extracting physical variables (scientific products) such as clock desynchronisation, electron 
content in the ionosphere (TEC), or range (instantaneous distance between the stations) from 
the code and phase measurements on ground and in space and auxiliary data (orbitography, 
internal delays, atmospheric parameters, ...). 

To this purpose we have developed the complete model of the time transfer at the re- 
quired 0.1 ps level. We have then developed in parallel two softwares: 
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1 . a program to simulate the raw MWL measurements which allows the addition of dif- 
ferent types of noise, biases, dead times in the measurements, phase ambiguities, etc 

2. an algorithm which implements the MWL model and provides the "scientific products" 
from the raw measurements. 

During the mission only the algorithm (2.) will be used, but the program (1.) is necessary 
for testing purposes. The two softwares are kept as independent as possible (different pro- 
gramming languages, different algorithms ...) to ensure a maximum efficiency of such tests. 
We provide some details on the software and the tests, considering different cases from the 
simplest to the more complex and realistic situation using real ISS orbitography data and 
MWL measurement noise from the MWL engineering model. 

The phase ambiguity removal of carrier phase measurements is performed by the algo- 
rithm and its success strongly depends on the noise of the observables. We have investigated 
the statistics of cycle slips which appear during this operation using experimental data ob- 
tained from the tests of the MWL engineering model. We present two novel methods which 
allow the reduction of the cycle slip probabilities by a factor greater than 5 compared to the 
standard method. 



1 Introduction 

Due to recent scientific breakthroughs such as laser cooling and atom trapping methods, 
huge progress has been achieved on the uncertainties of atomic clocks during the last twenty 
years (4j[8j ?,?]. Some of them reach a precision of a few parts over 10 17 in relative fre- 
quency. 

In a Space environment these atomic sensors would become an exceptional tool for 
promising applications in fundamental physics, geodesy, time and frequency metrology or 
in navigation. Onboard terrestrial or solar system satellites, their exceptional properties al- 
low them to test the fundamental laws of nature, to study the Earth's and solar system grav- 
itational potential and its evolution, or to explore the Universe 1211 . 

To this purpose, the ACES (Atomic Clock Ensemble in Space) mission 1151 . an ESA- 
CNES project will be installed onboard the ISS (International Space Station) in 2013. It 
consists of two atomic clocks and a two-way time transfer system (microwave link, MWL) 
with an overall uncertainty goal of 1 part in 10 16 after ten day integration (see section[2]for 
more details). This mission will carry out tests of fundamental physics such as testing the 
Einstein's Equivalence Principle or perform time transfer at the sub-picosecond level. 

This next generation of space clocks at the envisaged uncertainty level requires a fully 
relativistic analysis, not only of the clocks (in space and on the ground) but also of the 
time/frequency transfer method used to compare them [1 , 10, 13.20 5, 11). Similarly the 
modeling of the mission for the future data analysis requires a relativistic framework and 
the investigation of the complete model complying with the envisaged performances. 

The preparation of the mission data processing needs the development of two softwares. 
On one hand, because the mission is still in preparation, we need to simulate the MWL raw 
measurements in the most realistic way possible. On the other hand, the mission requires an 
algorithm to extract the "scientific products" from the raw measurements and the parameters 
of the mission. 

Finally, similarly to Global Positioning Systems, the phase ambiguity resolution is a key 
process which the mission success depends on. In fact an error of one cycle on one of the 
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three frequencies will lead to an error larger than the mission specifications. The occurrence 
probability of cycle slips strongly depends on the raw measurement noise and observable 
combinations must be found to reduce this probability. 

In this paper we study in more detail the preparation of the data analysis and provide 
some details on its testing which will allow to evaluate the high performances of the MWL. 
In sections [2] and [3] we briefly describe the ACES mission, and then the relativistic model 
used for the clocks and the time transfer to understand the nature of the raw measurements. 
We also relate their expression with all propagation effects between the two clocks (internal 
delays, atmospheric effects, relativistic contributions,...). This investigation is required to 
model the time transfer at the expected level. In section [4] we give information about the 
measurement simulation and the scientific products extracting algorithm and show some 
tests of the latter to evaluate its performances in terms of stability and accuracy. In the two 
last sections i[5]and[6]l before the conclusion, we discuss the phase ambiguity determination 
whose success strongly depends on the noise of the observables. We investigate the statistics 
of cycle slips which appear during this operation. To this purpose we use experimental data 
obtained from the tests of the MWL engineering model 1 16], and an ISS ephemeris. 



2 The ACES mission 

The ACES project led by the CNES and the ESA aims at setting up on the ISS several 
highly stable clocks around 2013. The ACES payload includes two clocks, a hydrogen maser 
(SHM developed by TEMEX) and a cold atom clock PHARAO (developed by CNES) re- 
spectively for short and long term performances, and a microwave link for communication 
and time/frequency comparison. The frequency stability of PHARAO onboard the ISS is 
expected to be better than 10~ 13 for one second, 3 • 10~ 16 over one day and 1 ■ 10~ 16 over 
ten days, with an accuracy goal of 1 ■ 10~ 16 in relative frequency. 
The ACES mission has as objectives : 

- to operate a cold atom clock in microgravity with a 100 mHz linewidth, 

- to compare the high performances of the two atomic clocks in space (PHARAO and 
SHM) and to obtain a stability of 3 ■ 10~ 16 over one day, 

- to perform time comparisons between the two space clocks and ground clocks, 

- to carry out tests of fundamental physics such as a gravitational redshift measurement 
and to search for a potential speed of light anisotropy and a possible drift of the fine 
structure constant a. 

- to perform precise measurements of the Total Electron Content (TEC) in the ionosphere, 
the tropospheric delay and the Newtonian potential. 

The time transfer is performed using a micro- wave two-way system, called Micro- Wave 
Link (MWL). An additional frequency is added in order to measure and correct the iono- 
spheric delay at the required level. It uses carriers of frequency 13.5, 14.7 et 2.25 GHz, 
modulated by pseudo random codes respectively at 10 8 s~ l , 10 8 s~ l and 10 6 s~ l chip rates. 
Moreover it has four channels that allow four ground stations to be compared with the ISS 
clock at the same time. 

According to the mission specifications, the microwave link has to synchronize two 
atomic clocks with a time stability of < 0.3 ps over 300 s, < 7 ps over one day, and < 23 ps 
over 10 days. The performance of this link is a key issue since it will perform high precision 
time comparisons without damaging the high performances of the clocks. 
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For our purposes we express the above requirements for the MWL in a simplified form 
by the temporal Allan deviation (o x (x)) : 

o 1 (t)=5.2-10" 12 -t"2 (1) 

for a single satellite pass over a ground station (for integration times x lower than 300 s) and 
by 

o x {t) = 2.4 -10" 14 -T3 (2) 

for longer integration times (for integration times X greater than 300 s). 

We take (fT]l and $2$ as our upper limits for all comparisons with the test results of the 
data processing algorithm in the following sections. 



3 The time transfer model 

Due to the outstanding performances of the clocks, some effects from fundamental physics 
have an impact on the clock behaviors and on the signal propagation. In this situation, a 
relativistic point of view must be considered to model the clocks and the time transfer. The 
aim of this section is to describe the relativistic model used for the ACES mission and 
the Micro-Wave Link. It is necessary to express the raw measurements as functions of the 
clock desynchronisation and of all effects which affect the signal propagation (relativity, 
atmosphere, internal delays, phase ambiguities, etc..) at the expected level. 

In a general relativistic framework each clock produces its own local proper time, in our 
case X s and X s for the ground and space clocks respectively. In order to model signal prop- 
agation between the ground and the space stations, we use a non-rotating geocentric space- 
time coordinate system. Thus t = xq/c is the geocentric coordinate time, x = (xi,X2,X3) 
are the spatial coordinates, where c is the speed of light in vacuum (c = 299792458 m.s ), 
We denote U(t, x ) as the total Newtonian potential at the coordinate time t and the posi- 
tion x with the convention that U > |19| . In these coordinates, the metric is given by an 
approximate solution of Einstein's equations valid for low velocity and potential << 1 

and £ « 1): 

ds 2 = _ (1 _ 2 J^ll )c 2 dt 2 + (1 + m±*l )<o* (3) 

where higher order terms can be neglected for our purposes 1201 . 

In this system, each emission or reception event (at the antenna phase center) is identified 
by its coordinate time t; (figure [5]l and a coordinate time interval is defined by Zy = tj — t-,. 
We define x g , ~v g and ~a g respectively as the position, the velocity and the acceleration 
of the ground station, and x s , v s and a s respectively as the position, the velocity and the 
acceleration of the space station. 

The /i frequency signal is emitted by the ground station at the coordinate time t\ and 
received by the space station at t^- The f% and frequency signals are emitted from the 
space station at tj, and fg , and received at the ground station at £4 and t(, . The third frequency is 
added to measure the TEC in the ionosphere which allows the correction of the ionospheric 
delay. 

According to equation (O, we introduce the notation [.] to express an interval of the 
coordinate time t in the proper time scale x A produced by a clock at x &{t). The following 
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Space station path 




Fig. 1 MWL principle 



equation relates a coordinate time interval T\2 elapsed between the coordinate time t\ and t 2 
to a proper time interval [T\i\ A elapsed in the proper time scale t a : 

where x a and va are respectively the position and the norm of the velocity of a point A 
expressed in the coordinate system described above, and the integral is evaluated along the 
path of the clock x 

The MWL is characterized by its continuous way of emission. It measures the time 
offsets between the locally generated signal and the received one. It provides three measure- 
ments (or observables) of the code (one on the space station, two on the ground) and three 
measurements of the phase of the carrier frequency at a sampling rate of one Hertz. 

An observable is related to the phase comparison between a signal derived from the 
local oscillator and the received signal, corrected for the frequency difference mainly due 
to the first order Doppler effect (see | 2 | for details of a similar procedure used in GPS). If 
we consider a special bit of the signal which is produced locally at x p and received at % a , an 
observable is a measurement of the local proper time interval between these two events. The 
observable is labeled with the arrival proper time x a . 

For instance, the upwards signal observable is 

4r ! «) = T^-< = T|-<. (5) 

The last identity comes from the fact that two identical segments of code are produced 
at the same local time : tJ = X s p . 
Thus we have 

AT S {T S (t 2 ))=T g ( tl )-T S (t 2 ). (6) 

We define the desynchronisation between the two clocks at the coordinate time t by the 
difference between the ground and the space proper times at t, T g (t) — T s (t). The raw mea- 
surement can be expressed as a function of the desynchronisation between the two clocks 
at the time of reception to, the emission and reception internal delays, and the propagation 
time between the two stations T\ 2 - For the sake of clarity and equation simplicity, we inten- 
tionally omit internal delays in the following expressions. They can be associated with the 
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propagation time except the fact they are proper time intervals in the local clock time scale, 
and their inclusion in the model poses no particular difficulties. 

For example, if we consider the observable from f\ signal PRN code, it is given by 

Ar s (x s (t 2 )) = - r ! (f 2 ) - [T n ]». (7) 

where T12 is the propagation time for the f\ signal. The latter is a coordinate time interval 
transformed to a ground proper time interval following equation ©. 

As we want to evaluate the desynchronisation between the two clocks, the time interval 
7j2 elapsed between emission from the ground station and reception by the satellite of the 
f\ frequency signal needs to be calculated. It can be written as 



_ gig | 2GM E ^ / x g (ti)+x s {t 2 ) +Rn 



where R\ 2 = | \Rn\\ = | \x~^(t 2 ) — x g (ti)\\, where the logarithmic term represents the Shapiro 
time delay 1171 (see e.g. (5) for a detailed derivation) and where A" 2 ° po and A'" 2 ° are respec- 
tively the tropospheric and ionospheric delays on the signal path. 

The phase observable for the f\ signal can be expressed similarly except for two impor- 
tant features : the effect of ionosphere and the phase ambiguity. The ionosphere delay takes 
opposite signs for code and phase [3|. Moreover a carrier phase measurement is less noisy 
than code measurement but affected by a phase ambiguity. This ambiguity corresponds to 
integer number Nj of the signal period 1 /fi. For instance the phase measurement of the f\ 
signal is given by 

at; (T s (t 2 ))= T «( f2 ) - T s (r 2 ) - pr£]« (9) 

where Tf 2 ls the phase propagation time of the f\ signal and given by ® with a change 
of sign of the ionospheric term A'\ 2 no ■ 

An error in the determination of is referred to as a "cycle slip". It will be studied in 
detail in sections <[5j> and ©. 

The expressions of code and phase observables for the two remaining frequencies can 
be derived similarly to I0 and The combination of f\ and f 2 signal observables gives 
the expression of the desynchronisation between ground and space clocks : 



THt a )-x s {t a ) = -\Ax\x\t 2 )) -4T«(T«(r 4 )) 

+ T\2 — ?34 

_ p- { u( t ,r s ) , v 2 g ( t ) 



(10) 
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where t a = '■ 2 ^ A , and where Ax s (i s (t2)) and Ax 8 (z g (t4j) are the observables respec- 
tively from the ground and onboard the satellite, and where we have neglected non-linearities 
of T s (t) and T s (r) over the interval £4 — (few milliseconds). The integral terms result from 
proper time to coordinate time transformations. They are small corrections of order 10~ 12 s 
to the desynchronisation I s {t a ) — X s (t a ). 

In d 10b . the difference T\2 — T14 needs to be calculated from the knowledge of satellite 
and ground positions and velocities obtained from orbit restitution (see |6j). It includes the 
determination of ionospheric delays through the evaluation of the Total Electron Content 
(TEC). This is done by combining the observables coming from fx and fa signals because 
the signal paths are approximatively the same and there is almost one order of magnitude 
between the two frequencies which provides a good estimation of the TEC. 

The same equation exists for phase observables. Considering phase ambiguities are cor- 
rectly removed, it differs from the previous equation d 1 fr by the term T\i — T34. This term 
is not the same for code and phase desynchronisation equation due to the different effects 
from ionosphere on signal propagation. 

4 Data analysis and simulation 

The data analysis of ACES measurements has specific goals. It aims at providing the searched 
scientific products from the MWL measurements and some additional inputs, in particular 
ISS orbitography, parameters for the troposphere model, etc... Its goal is to evaluate physi- 
cal parameters which will serve to study the Earth's atmosphere or to probe the fundamental 
laws of physics. 

However, the data analysis has to deal with experimental issues such as noise on mea- 
surements or loss of signal acquisition which leads to dead times in the observables. Even 
in adverse conditions, it has to provide results within the expected performances of the mis- 
sion. Thus the data analysis must be prepared before the launch of the mission, and tested in 
all pessimistic situations. 

To this purpose, two software packages have been developed considering the complete 
model of time transfer described in (6J. On one hand, a raw MWL measurement simulation 
which produces the observables from given clock difference, orbit, internal delays, iono- 
sphere and troposphere parameters, and on the other hand, a data analysis algorithm which 
extracts the searched parameters from the raw measurements. The two softwares have been 
implemented in the aim of keeping them as independent as possible. That is why they use 
different languages of programmation and different algorithms of calculation. 

In the following paragraphs we will give some details on both softwares. 

4. 1 The raw measurement simulation 

The simulation produces the raw measurements of the MWL for code and carrier phase on 
the three frequencies. To this purpose it needs several parameters chosen by the user such 
as the clock behaviors (T s (t) and T ! (f)) , the space and ground station trajectories or the 
evolution of internal delays. 

It is implemented in an object oriented language and simulates the six observables using 
an iteration procedure. 

As an example, the simulation of the f\ signal code measurement At s (z s fa)) measured 
at the space proper time T s (t2) starts with the calculation of the propagation time T between 
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the ground and space stations respectively located at x g (f2) and x s (fc), which are obtained 
by interpolation of the orbitography data. Due to the propagation delay, the ground station 
was not at x g (t2), but at a previous position along its trajectory when emitting the signal. 
New position, velocity and acceleration are interpolated at the coordinate time t2 — T which 
allows to calculate a new value for the propagation time T. By convergence, this method 
provides the true time of signal emission t\ by alternatively calculating the propagation time 
T for a signal arriving at the coordinate time t<i and interpolating the ground station position, 
velocity and acceleration. We then evaluate the proper time T g (?i) of signal emission from 
the ground station and finally the code observable 4t s (t" 5 (?2))= T g (h) ~ t s {h)- 

Once the "theoretical" measurement is calculated, perturbations are added to it. The 
simulation allows addition of measurement noise, clock error, internal delays and errors on 
them, antenna phase patterns, dead times in measurements, orbitography error, etc ... 

4.2 The data processing algorithm 

The MWL data analysis implements the MWL model and provides the scientific products 
from the raw MWL measurements. Its outputs are the Total Electron Content (TEC) in 
the ionosphere along the line of sight, the desynchronisation between the ground and the 
space clocks, the tropospheric delay and the instantaneous distance between the two stations. 
These products are obtained from code only, and code plus phase measurements. Besides 
the raw measurements, the calculation needs some other data such as the orbitography of the 
ground and space stations, atmospheric parameters (temperature, pressure, ...) or the internal 
delays. 

Contrary to the raw measurement simulation, the data processing algorithm is written in 
an assembly language. Before the calculation of the scientific products, it performs a prepro- 
cessing so as to detect dead times in code and phase measurements, to combine correctly the 
measurements of the three frequencies and to transform space and ground station trajecto- 
ries in the adapted coordinate system. Moreover it also determines cycle ambiguities on the 
carrier phase measurements before using the non-ambiguous phase observables to estimate 
the scientific products. 

The calculation is based on analytical expressions at a chosen time. For instance, in 
expression fllOl l, the term T\% — T34 is evaluated with a Taylor expansion. Its expression is a 
function of position, velocity and acceleration of the two stations at one chosen coordinate 
time. The same kind of expansions are used for the calculation of TEC, range or tropospheric 
delay. 

4.3 Tests and results 

The testing of the algorithm is performed by taking into account independently and then 
simultaneously the different noises and perturbations which will appear during the mission 
(code and phase measurements noise, clock noise, orbit restitution noise, dead times, ...). It 
allows to verify its stability and accuracy performances. 

To this purpose we feed the output of the raw measurement simulation into the data 
analysis algorithm and compare if its output corresponds to the initial given functions used 
as input for the simulation. The differences between the outputs of the algorithm and the 
initial functions are identified as final errors on scientific products. In the actual mission, 
only the MWL data analysis will be used, to obtain directly the required scientific products. 
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Fig. 2 Error on desynchronisation when no noise is added 
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Fig. 3 Temporal Allan deviation when no noise is added 



Many tests have been carried out, but for clarity, we only show the results of three of 
them. First we add no noise which would damage the performances of the link. The ob- 
tained error on desynchronisation is drawn on figure [2] and stays under 0.1 picoseconds. 
The remaining term corresponds to difference of tropospheric and ionospheric delays for 
Ku signals (f\ and /a) which have not crossed the same atmospheric layers. The resulting 
deviation of the error on the desynchronisation is two orders of magnitude below the spec- 
ifications (cf. figure [3}. We conclude that no term giving a deviation over the specifications 
has been neglected. 

Now the three observables have dead time intervals. They can be caused for example 
by a loss of station visibility or to the occultation of the reception antenna by a solar panel. 
Figure [4] shows that the software deals with dead times and provides the scientific prod- 
ucts without error. The corresponding deviation is also two orders of magnitude under the 
specifications of the mission. 

Finally the last presented test is performed adding white noise on all simulated raw mea- 
surements with an amplitude corresponding to the specifications (l}. The resulting deviation 
stays under the specifications (cf. figure[5]l. It is calculated in accordance with equation < | 1 b 
: the white noises from f\ and f% observables are not correlated which explains the square 
root of two gain with respect to the mission specifications. The addition of noise is han- 
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Fig. 4 Error on desynchronisation with dead times in measurements 
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Fig. 5 Temporal Allan deviation considering white measurement noise on observables 



died by the data analysis software and this is demonstrating the robustness of the algorithm 
confronting this difficulty. 

As written before, many other successful tests have been done considering different 
perturbations affecting the accuracy and the stability of the ACES time transfer. They involve 
realistic orbitography error, random walk noise on clocks or offset and noise on internal 
delays. The algorithm survived all these tests which complies with its initial objectives. In 
addition, the combination of the two softwares allows one to easily evaluate the effects of 
each parameter on the time and frequency transfer or to find out what the final results would 
be if there was an Einstein Equivalence Principe violation. 



5 Phase ambiguity resolution 

For all time/frequency transfer techniques the carrier phase observables have the advantage 
of significantly better precision than the code measurements, and therefore their use has 
been implemented in most techniques. Using the phase observables requires resolution of 
the phase ambiguities, which has been investigated in detail in many techniques, in particular 
the Global Positioning System |7]|9). 
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However there are several differences between GPS and the Micro Wave Link. On one 
hand, the MWL is a two way time transfer with three frequencies involved. The combina- 
tion of measurements coming from an upwards and a downwards signal allows to eliminate 
restraining factors without canceling the clock behaviors. On the other hand, the frequen- 
cies used in the MWL are higher than GPS frequencies except for the third MWL frequency 
which is added to evaluate the TEC in the ionosphere. There is almost one order of mag- 
nitude between the lowest and the highest frequencies. The phase ambiguity identification 
must therefore be investigated for the particular case of the ACES mission. 

Similarly to equation ( |10| l, the desynchronsation between the two clocks can be ex- 
pressed as a function of phase observables and cycle ambiguities. It allows to study the 
consequences of the phase ambiguity detection on the final calculation. 



where the last term is the correction for the leading order ionospheric effect and A 
includes terms due to path asymmetries, tropospheric delays, relativistic corrections, higher 
order ionospheric terms, etc... (see (6) for details). Although these terms need to be taken 
into account, they are not relevant for the discussion concerning phase ambiguity resolution 
and will be neglected in the following. 

According to equation d 1 1 b - the cycle slip detection is necessary for the three frequen- 
cies. In fact an error of one cycle on f\ , fi and (i.e. an error of one unit in the deter- 
mination of N\, N2 and N3) leads respectively to an error of 3.7- lfT 11 s, 3.4- 10" 11 s and 
1 .0 • 10~ 12 s on the clock desynchronisation, significantly larger than the uncertainties aimed 
at. 

The carrier phase ambiguity Nj on the /, phase observable is evaluated by combining the 
phase and code measurements and fixing Ni to the nearest integer. These unknown numbers 
of periods added to phase measurements stay constant for continuous data ie. as long as the 
signal acquisition is kept. This means we obtain during a pass several evaluations of the same 
Ni and then we can average the noise on them to get a more accurate value. Consequently 
correct cycle slip identification depends strongly on the observable noise types and levels. 

In order to express the phase ambiguities as a combination of observables, we have to 
use a number of hypothesis. First the path of code and phase measurements tagged with the 
same coordinate time and at the same frequency are identical. We neglect the change on 
geometric delays due to the opposite effect of ionosphere on code and phase propagations. 
Moreover we assume all downwards signals received at the same coordinate time by the 
ground station have the same path through the atmosphere. Finally the part of ionosphere 
crossed by upwards and downwards signals are supposed identical for combined observables 
(i.e. for signals within a few ms of each other). 

When forming the difference between phase and code observables measured at the same 
coordinate time and at the same frequency, we cancel all the propagation terms except the 
ionospheric delay (see equations ©, ^ and $9$). Actually the ionospheric term depends 



r g (t a ) - t s (t a ) = \Ut; (x s {t 2 ))-Axl (T«(r 4 )) 



fi 




(11) 
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on the signal frequency and on the measurement nature, code or carrier phase. The leading 
term is proportional to the inverse of the square of the considered frequency. In addition, for 
the same frequency, the effect of ionosphere on code is at first order the opposite effect of 
ionosphere on carrier phase. The error on the ionospheric delay comes from the error on the 
TEC determination and is dominated by its leading term. Then the combination of equations 
|7ll and $9$ for the ambiguity Ni on the first frequency can be written as : 

^=Ar;(r\t 2 ))-AT°(r\t 2 ))-2-^^TEC. (12) 

The same kind of equation can be written for f 2 using code and phase observables, 
A T 8 ( T g (f4 ) ) and A t| ( T g (?4 ) ) , or for the third frequency fo with A T g ( T g {t& ) ) and A t| ( T g (t$ ) ) . 
They relate the phase ambiguity N with the difference between phase and code observables 
and twice the ionospheric delay at the considered frequency. Consequently the error on the 
phase ambiguity 8N, depends on the code and phase measurement noises, respectively 5C,- 
and 8<Pi, and on the error on the TEC evaluation. The latter is determined with the code 
measurements from f 2 and f$ signals [6]: 

^^TEC=-^^AT g (T g (t 4 ))-AT g (T g (t 6 ))"j, (13) 

The Total Electron Content is evaluated with this formula and the error on its deter- 
mination is dominated by f 2 and fj code measurement error. Hence we obtain, for each 
frequency fi, the dependence of the error on the phase ambiguity SNi on the uncertainty in 
the observables by combining equations | |12| | and d 1 3b : 

for/i : ^- = 54>i-5Ci-5.7-10" 2 (5C 2 -5C3), (14) 
fi 

for f 2 -^P L = 8<P 2 - 8C 2 -4.8 • 10" 2 (5C 2 - 8C3), (15) 
12 

8N3 

for fi : — - = 54> 3 - 8C3 - 2.04 • {8C 2 - 8C3). (16) 

J3 

In the last equation d 1 6b . two terms involve the fj code noise 8C3. They are perfectly 
correlated and partially cancel each other. 

The noise model is evaluated with the latest results from the engineering model [16|. 
These results correspond to code and phase measurements performed on a f\ signal for a 
received power of -95 dBm (this S/N ratio is related to a zenithal position of the ISS over the 
ground station at about 350 km). They show respectively a standard deviation of 1 • 10~ 12 s 
and 1 ■ 10~ 13 s for the code and the phase measurements. This 1 s noise will be used in this 
section to estimate the statistics of cycle slips. 

For this purpose, we assume the noise on f 2 approximatively the same. On the contrary, 
the noise on fa is assumed to be 7 (~ /1//3) times larger on the carrier phase and 100 
(ratio of chip rates) times larger on the code. In the following table (Tab. [3, we give the 2a 
uncertainty (at 1 s points with no averaging) of each term appearing in the equations d 1 4b - 
( 115b and d 1 6b . These noise levels have to be compared with the half of the period 1 / fi due 
to the fixing of the ambiguity Ni to the nearest integer : if the 2a noise level is equal to half 
the period, there is a 5 % chance of a cycle slip. 

We note that at all frequencies the cycle slip probability is less than 5 % - remember that 
on /3 the code noise and half the ionosphere noise cancel (see equation i ll 6b ) - even when 
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Table 1 95 % confidence (2o~ uncertainty) of terms appearing in the phase ambiguity resolution according to 
equations )14t . )15t and i\6\ for an input power of -95 dBm 



95 % confidence (2o~ uncertainty) 


/l 


h 


h 


Noise on code (s) 
Noise on carrier phase (s) 
Ionospheric error due to the noise on code (s) 
Signal period \ j f\ (s) 


2-lfT 17 
2-l(T 13 
1.1-10" 11 
7.4- 10-" 


2-lfT 17 
2-l(T 13 

9.6- 1(T 12 

6.7- 10"" 


2 .10- 1U 
1.4- IQ- 12 
2-2-KT 10 
4.5 10" 10 



using individual 1 s measurements (no averaging). The probability can be further reduced 
by averaging the measurements of a continuous observation. 

However two facts will disturb the cycle slip identification. First the received power does 
not remain at -95 dBm but decreases when the distance between the two stations increases. 
As an example, when the space station is at a ten degrees elevation, the received signal 
power is -115 dBm. From the tests of the engineering model, this corresponds to an eight 
fold increase of the code noise. In these worst conditions, the direct method of ambiguity 
resolution fails. One solution is to carry out a weighted average with all 1 s measurements 
of a given continuous pass. Indeed, the noise level of each measurement will depend on the 
incoming power (varying as a function of elevation). The dependence of the noise on power 
has been measured so the weighting can be preformed as a function of the locally measured 
power level. 

Secondly the measurement noise is not perfectly white. Some flicker noise damages the 
averaging and the ambiguity removal, and even averaging over the complete pass is insuffi- 
cient due to the presence of flicker noise. We need to search for observable combinations to 
increase the ambiguity resolution success rate. 

The figures of Tab. [T] show that the leading error terms come from the code noise on the 
third frequency. Our aim is to get rid of the code measurement of fa, or at least to reduce its 
impact on the final statistics. 

We start by noting that, according to the noise levels in Tab.[T] the ambiguity resolution 
for the Ku signals {f\ and fa) is more likely to be successful than for the third frequency. 
Then, if we suppose these ambiguity resolutions are correctly performed (this is likely the 
case, see section[6]l, these considerations bring two new methods for the fa signal ambiguity 
resolution. On the one hand, the combination of fa and fa phase measurements gives a new 
expression of N$ which does not depend on the fa code measurement : 



= A^(t 6 ))-A^(t 4 )) 



N 2 fj-fj 40.308 

+ h + m c TEC - 



(17) 



In the previous equation, the Total Electron Content is evaluated by the combination of 
fa phase and code measurements : 

^rac=|L(t*w)-4T»(t»« ( ))-|j. (18) 

It allows a new evaluation of the TEC whose error depends only on 8C2 and 8<p2 and 
which can replace the equation d 1 3 b - Inserting equation i ll 8[ i into H7[ leads to 
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— ± = 803 -22.72 -5*2 +21.72 SC 2 . (19) 

J3 

On the other hand, the second method is to substitute d 1 3b into ( 112b (but written for fa) 
for half the ionospheric term and using dl8l l for the other half. The error from the fa code 
measurement almost cancels, and the error on N3 is : 

O IT 

— = 54>3+2.3-10" 2 -5C 1 -22.22-54>2+21.20-5C 2 . (20) 

J3 

Equations dl9b and d20b reduce the error from the code measurements from fa and then 
increase the success rate. The choice between one of these methods depends on the real 
mission noise levels. With the noise measured in the engineering model tests, they both 
approximatively give a five fold gain on the uncertainty of 8N3 . 

In conclusion, the statistics of successful carrier phase ambiguity resolution strongly de- 
pend on the pass characteristics and its relation with the noise levels. These kind of realistic 
received noises are studied in the following section[6] 



6 Tests with expected noise 

In this section, we investigate the error on the phase ambiguity resolution from realistic 
measurement noises and try to reduce the probability of cycle slips on the ambiguity iden- 
tification. For an ISS pass over a ground station -ie. for time duration less than 600 S-, 
the noise measured on the Engineering model can be modelled in a simplified from by its 
temporal Allan deviation. We use respectively for code and phase measurements at -95 dBm 
(o c x {x)f= (1 • 10-%-:? ) 2 + (2 • 10- 13 ) 2 s and (af (t)) 2 = (1 • 10-%-? ) 2 + (7 • lO" 14 ) 2 s, 
where % is the integration time [ 16 1. 

The measured noise is composed of white noise and flicker noise. The former is averaged 
rapidly as the square root of time. Then the noise average is limited by the flicker noise. On 
figure ©, we compare the averaging of the measured noise on fa code at -95 dBm, with the 
averaging of a pure white noise with the same level at one second (a = 1 • 10~ 12 s). 

During a pass of the International Space Station over a ground station, the signal over 
noise ratio strongly varies and its received signal power can decrease by over 20 dBm. This 
ratio actually depends on the distance between the two stations, the elevation of ISS and 
others parameters (antenna orientation, weather conditions, etc ...). A model provided by 
EADS-Astrium gives power evolution with distance and elevation. The engineering model 
measurements of the noise dependence on power carried out at TimeTech 1 16] then allow 
us to derive the dependence of the noise on elevation and distance, i.e. the variation of the 
noise during a given pass. 

We perform Monte Carlo simulations to estimate the probabilities of a cycle slip and 
its uncertainties (95 % confidence). In fact, we create code and phase signals with the same 
noise properties as the measured noise on the engineering model for the frequency fa . The 
fa and fa noise for code and phase measurements are similarly generated taking into ac- 
count chip rate or frequency ratios. Then these noises are modulated as a function of the 
elevation and distance for a given pass. In a pessimistic approach we suppose the two kinds 
of noise (white and flicker) on code and phase measurements are affected in the same way 
by the received power. Hence we can calculate the 95 % confidence of the ambiguity error 
in accordance with equations d 14b . dl5b . dl6b . l |19b and ( 120b . 
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Fig. 6 95 % confidence of averaging of white noise (red) and f\ code measured noise at -95 dBm (blue). 
Both have the same deviation at one second. 

The noise on the ambiguity identification is expressed in seconds : as this noise has to 
be compared with the signal period, we express it in term of a fraction 8Ni of the consid- 
ered signal cycle. We average over a number of 1 s points to decrease the uncertainty on 
Nj. Then, if 8Nj after averaging exceeds half a cycle with more than 5 % probability, we 
consider that the ambiguity resolution has failed as N, is fixed to the wrong number. As 
an improvement over simple averaging, we also consider weighted averaging, where the 
weights are determined by the noise level of each point (depending on received power, ie. 
elevation and distance). 

In the following we study the time transfer between the International Space Station and 
a ground station based in Toulouse, France (43°36'N, l°26'E ) to simulate realistic noise lev- 
els. Actually, this station has been chosen as the master ground station of the ACES mission. 
Then we choose an ISS trajectory corresponding to one pass over the ground station. To this 
purpose, we consider an ephemeris of ISS corresponding to the 20 ,h of May, 2005. How- 
ever in order to simulate several passes with different elevation higher than four degrees, we 
slightly shift the time origin of the ground station trajectory. This way we obtain numerous 
passes characterized by their reached maximum elevation and whose duration is between 
350 and 600 s. 

Figure l|7]l shows the uncertainties (95 % confidence) of 8N\ and 8N2 as a function of 
the maximal elevation reached by the space station during the pass. These uncertainties are 
calculated with and without a weighted averaging over a complete pass ie. data are acquired 
continuously during the pass. It allows to assess the effect of weighting on N\ and N2 final 
estimation. 

Weighted averaging leads to a gain on the success rate and particularly with high eleva- 
tions : at ninety degrees elevation, the statistics are almost divided by a factor six between 
the two averaging methods. In fact, at lower elevations, the noise level on the observables 
does not sensibly change during the pass, and in this case the weighting becomes ineffi- 
cient. We note that above six degrees maximum elevation, the phase ambiguity is correctly 
determined (95 % confidence). 

Finally we investigate the impact of averaging (arithmetic or weighted) and calculation 
(in accordance with equations d 1 6b . il9l or d20l l) methods on the uncertainty on the phase 
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Fig. 7 95 % confidence of the phase ambiguity resolution error as a function of the maximal elevation of the 
pass. The four curves show this dependence for f\ (blue) and fa (red) signals using arithmetic (cross-dashed) 
or weighted (solid) averaging in accordance with equations )14t and )15t . 

ambiguity N3. For the noise considered in this paper, the methods coming from equations 
l !19t and < |20b give approximatively identical results. 
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Fig. 8 95 % confidence of the phase ambiguity resolution error for the fj signal as a function of the maximal 
elevation of the pass. The four curves show this evolution whether we consider first equation )16t (cross- 
dashed) or proposed method ({19} or {20} ) (solid) and whether arithmetic (red) or weighted (blue) averaging 
is performed. 

Figure [8] depicts the different combination of approaches concerning the averaging and 
the calculation. Contrary to weighting, the observable combinations given by equation d 1 9b 
or l !20| > brings a factor five gain for all elevations (see the two lowest curves on figure[8]>. We 
also remark that the proposed measurement combinations H9i or d20b are essential in the 
cycle slip detection on the third frequency, and thus for the mission objectives. 
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However the obtained probability of successful ambiguity resolution on needs correct 
identification of the Ku cycle slips {f\ and fj). Consequently, an unbiased estimation of the 
scientific products is not limited anymore by the phase ambiguity measurement but by 
the two first frequencies. But in any case, using weighted averaging and methods dl 9b and 
d20b for /3, the cycle ambiguity resolution has less than 1 % probability to fail for any 
uninterrupted satellite pass with a maximum elevation over four degrees if we suppose that 
fi and fi phase ambiguities have been correctly solved. 

7 Discussion and conclusion 

The models of the time transfer of the ACES mission and the associated relativistic effects 
were investigated at the 0. 1 picosecond level. In addition, the measurements were related to 
all physical effects appearing during the station communication. 

Then these models were applied to develop an end-to-end simulation. In fact, two soft- 
wares were implemented : a simulation which creates the raw measurements of the mission 
and allows addition of different kinds of perturbations, and an algorithm which uses these 
raw measurements to extract scientific products. 

Furthermore the MWL phase ambiguity resolution was tackled and methods were stud- 
ied to reduce the failure rate of this process. The cycle slip removal faced two difficulties, 
the presence of Flicker noise which can not be averaged, and the evolution of the signal over 
noise ratio during the pass. Considering a model for the power evolution with distance and 
elevation, we simulated noises complying with experimental code and phase measurements 
on the first frequency f\ . Then the fi and fa code and phase noises were deduced assum- 
ing they are related to fi code and phase noise with respectively chip rate and frequency 
ratios. In a pessimistic approach, we supposed all components of noise (white and Flicker) 
are similarly modulated by the S/N ratio. 

With Monte Carlo simulations, we obtained the most realistic estimates of the result- 
ing performance for uninterrupted passes. The success rate is now limited by the ambiguity 
resolution of the Ku-band signals. In fact the high probability of successful resolution on 
the third frequency is based on the hypothesis that there is no cycle slip on f\ and fa- Their 
resolutions are limited by the code noise from frequency fa . The failed resolution proportion 
for Ku-band signals is less than 5 % for passes with maximum elevation over six degrees. 
Consequently, in this case, the phase ambiguity resolution is successful for all three fre- 
quencies with more than 94 % probability. An increase of the S-band signal chip rate, from 
10 6 i _1 to 5 • 10 6 i _1 , would strongly improve the Ku-band signal phase ambiguity as it will 
approximatively divide fa code noise by a factor 5. 

On one hand, these statistics are based on a pessimistic hypothesis of noise variation 
with received power. We assume indeed all kind of noise in the measurement are similarly 
modulated by the power evolution whereas it is likely that the Flicker noise is not or less af- 
fected. On the other hand, we have optimistically considered only completed uninterrupted 
passes over the ground station. In a close future, we will investigate the phase ambiguity res- 
olution for incomplete passes (loss of lock during the pass) and other perturbations affecting 
the raw measurements during a pass (eg. temperature variations). 
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